AD-A168  362 
UNCLASSIFIED 


METRIC  CALIBRATION  OF  THE  MILLSTONE  HILL  L-BANO  RADAR 
<U>  MASSACHUSETTS  INST  OF  TECH  LEXINGTON  LINCOLN  LAB 
E  M  GAPOSCHKIN  19  AUG  85  TR-721  ESD-TR-85-172 
F19628-85-C-B8B2  F/G  1779 


1/1 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  RDNCAu  Of  STANOANOR- l»»J- A 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 


METRIC  CALIBRATION 
OF  THE  MILLSTONE  HILL  L-BAND  RADAR 


E  M.  GAPOSCHK1N 
Group  9J 


TECHNICAL  REPORT  721 
19  AUG1ST  1985 


DT1C 

SELECTEt|\ 

OCT  1 7  ©85  j  I 

B 


Approved  Tor  public  release;  distribution  unlimited. 


LEXINGTON 


MASSACHUSETTS 


Abstract 


The  precision  Lageos  orbit  is  used  for  near  real  time 
calibration  of  the  Millstone  Hill  Radar  1MHz  chirped  frequency 
domain  tracker.  With  independent  tracking  data,  the  range 
calibration  is  given  to  «  lm,  with  Millstone  data  the  range 
calibration  is  given  to  «  5m.  The  range  rate  is  an  unbiased 
measurement,  with  an  accuracy  of  «  liran/sec .  The  azimuth  and 
elevation  are  calibrated  to  0°.020  and  0#.002  respectively. 
Tables  of  calibration  constants  are  given. 
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1.  Introduction 

Calibration  of  metric  data  from  the  Millstone  Hill  Radar 
requires  an  external  standard.  The  basic  observables  (time 
delay,  delay  rate,  etc)  must  be  converted  into  metric  quantities 
using  a  calibration  model.  The  parameterization  of  the 
calibration  model  (i.e.  the  transformation  from  observables  to 
metric  quantities)  comes  from  the  system  designer.  The  numerical 
values  used  in  the  calibration  model  come  from  the  calibration 
against  an  external  standard.  We  review  here  the  calibration 
model,  data  reduction  methods  used  at  the  Millstone  Hill  Radar, 
and  the  methods  used  to  obtain  the  calibration  coefficients. 
Quantities  that  are  measured  within  the  system,  and  used  to 
determine  the  observable,  such  as  peak  power,  bandwidth,  and 
droop,  are  not  considered  part  of  calibration.  The  external 
standard  is  a  satellite  (Lageos,  1976  039  01,  SDCI8820)  with  a 
very  well  known  orbit.  Therefore  the  methods  used  in  this 
analysis  can  be  applied  to  other  sensors.  The  limitations  of 
this  approach  will  also  be  described.  Because  the  metric 
accuracy  of  the  radar  is  continually  being  improved,  so  too  are 
the  calibration  techniques.  However,  we  expect  this  method  to 
provide  the  most  accurate  independent  calibration  of  a  metric 
radar.  We  describe  here  the  present  status  and  results  from  this 
calibration  effort,  in  preparation  for  the  next  overall 
improvement  in  system  performance. 


2.  Objectives 

The  Millstone  Hill  Radar  is  used  primarily  for  space 
surveillance  of  objects  in  orbit  around  the  earth.  The  radar 
operates  at  a  fundamental  frequency  of  1295  MHz,  and  can  observe 
targets  from  close  earth  orbit  (several  hundred  kilometers  range) 
to  geosynchronous  orbit  (6  earth  radii  or  42000  Km)  and  beyond, 
depending,  of  course,  on  radar  cross  section.  Recent 
modernization  and  refined  signal  processing  methods  (1)  give  an 
expected  precision  of  better  than  lm  in  range  and  lmm/sec  in 
range  rate,  for  Lageos.  The  present  mechanical  azimuth  and 
elevation  encoders  can  provide  an  angle  precision  of  0°.002.  The 
task  of  calibration  then,  is  to  transform  this  precision  into 
accuracy. 

In  addition  to  the  need  for  high  accuracy,  there  is 
timeliness.  Data  taken  with  the  radar  is  generally  used  within 
minutes.  Therefore,  the  calibration  would  best  be  available 
immediately.  Of  course,  refined  calibration  constants  are  always 
useful  in  a  posteriori  analysis. 

The  reduction  of  basic  observables  (e.g.  time  delay)  to 
metric  quantities  (e.g.  range)  also  requires  that  certain 
environmental  corrections  be  made.  The  most  important  is 
tropospheric  and  ionospheric  refraction.  This  study  describes 
the  simple  models  used.  It  is  quite  likely  that  limitations  in 
calculating  refraction  imposes  the  ultimate  limit  on  absolute 
accuracy. 
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In  summary  we  intend  to  provide  near  real  time  calibration 
of  the  Millstone  Hill  Radar  to  one  meter  and  one  millimeter  per 
second,  or  better. 

In  addition  to  providing  absolute  calibrated  metric  data,  in 
near  real  time,  the  calibration  procedure  also  provides  the 
following  capabilities: 

1)  Monitor  of  system  performance. 

2)  Identify  and  flag  error  conditions,  e.g.  clock 
malfunction . 

3)  Monitor  and  assess  environmental  parameters. 

4)  Provide  rapid  assessment  of  system  modification, 
enhancements,  etc. 
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3.  Method 


Our  independent  yardstick  for  calibration  is  the  Lageos 
satellite.  The  orbit  of  Lageos  (Fig.  1)  is  the  best  known  orbit 
in  the  solar  system,  that  is  to  say  it  can  be  known  in  a 
geocentric  inertial  reference  frame  with  an  accuracy  of  better 
than  20cm.  The  orbital  and  physical  characteristics  of  Lageos 
are  given  in  Table  I. 


TABLE  I 

LAGEOS  CHARACTERISTICS 


Semi-major  axis 

eccentricity 

inclination 

Mean  Motion 

Shape 

Mass 

radius 

surface 

core 

reflectivity  coeff. 
laser  ranging 

radar  cross  section  (L-Band) 


12000km 
0.006 
109°. 0 
6.6  rev/day 
sphere 
407.00  kg 
30.0  cm 
aluminum 
brass-copper 
1 .002 

426  cube  corner  reflectors 
0.28  m2 


EQUATOR 


These  characteristics  make  Lageos  an  excellent  calibration 
target  for  a  number  of  reasons: 

1)  The  altitude  and  eccentricity  reduce  the  effects  of  the 
anomalous  geopotential.  The  small  area  to  mass  ratio 
reduces  the  effects  of  non-gravitational  forces.  The 
spherical  shape  makes  these  forces  easier  to  compute. 
From  an  orbital  accuracy  standpoint,  it  is  optimized. 

2)  Lageos  has  a  good  radar  cross  section  at  L-Band,  and 
precision  metric  observations  are  possible. 

3)  Lageos  has  good  station  geometry,  which  allows 
determination  of  azimuth  and  elevation  dependent 
calibration  parameters. 

4)  The  laser  cube  corners  allow  observations  with  laser 
ranging  systems.  There  is  a  global  network  of  laser 
tracking  systems  with  tracking  accuracy  of  better  than 
10cm.  These  observations  allow  determination  of  the 
Lageos  orbit  to  comparable  accuracy. 

While  the  laser  ranging  systems  have  the  capability  to 
provide  data  with  an  accuracy  of  better  than  10  cm,  it  requires  a 
number  of  post  processing  corrections  to  achieve  this  accuracy. 
Laser  ranging  data  is  available  within  minutes  of  the  observation 
without  these  refinements,  and  is  known  as  quick  look  data.  The 
accuracy  of  quick  look  laser  data  ranges  from  10  cm  to  50  cm.  It 
is  the  quick  look  data  that  we  exploit  in  the  real  time 
calibration  system. 


There  are  four  modes  of  using  Lageos  for  calibration: 

1)  Based  on  the  metric  data  alone,  an  orbit  is  determined 
with  calibration  parameters,  i.e.  simultaneous  recovery  of  the 
state  vector  at  some  epoch,  and  the  sought  calibration 
parameters.  This  is  rather  time  consuming  (slow)  in  that  it  will 
use  data  from  many  passes  over  several  days  to  get  the  benefit  of 
orbital  constraints,  and  a  complete  force  model  is  needed  to 
achieve  submeter  accuracy.  Therefore  it  is  not  suitable  for  near 
real  time  application. 

2)  With  a  state  vector  derived,  say  as  in  1)  above,  the 
orbit  can  be  predicted  to  the  time  of  observation.  This 
prediction  can  be  compared  with  the  actual  observation,  and  the 
difference  can  be  used  to  obtain  calibration  parameters.  This 
method  is  limited  by  the  degradation  of  orbital  accuracy  as  the 
observation  time  gets  far  from  the  epoch  or  data  interval  used  to 
obtain  the  state  vector.  However  it  is  fast. 

3)  Based  on  independent  tracking  data,  say  laser  ranging, 
a  precision  orbit  is  determined.  Then,  the  differences  between 
that  reference  orbit,  and  the  observation  can  be  used  to  obtain 
calibration  values.  This  approach  involves  heavy  computation  as 
in  1)  above.  It  has  the  advantage  of  a  more  accurate  orbit,  and 
independent  determination  of  orbit  and  calibration  parameters. 
This  is  the  most  accurate  method. 


4)  With  a  state  vector  derived  in  3)  above,  one  can 
proceed  as  in  2).  This  has  the  same  relative  advantages  and 
drawbacks.  The  orbital  accuracy  is  better,  and  the  extrapolation 
accuracy  will  also  be  better. 

The  basic  results  reported  here  are  derived  from  1,  2,  and  3 
above.  Mode  2  (and  Mode  4)  have  been  implemented  in  a  software 
package  described  by  Kantrowitz  (2)  and  Byers  (3)  and  a  few 
remarks  about  it  are  in  order.  The  precision  orbit  computation 
is  described  in  section  6  below.  The  program  (DYNAMO)  has  an 
output  option  to  generate  predicted  observations  (azimuth, 
elevation,  range,  range  rate)  from  an  initial  state  vector.  The 
predictions  are  made  for  every  pass,  typically  for  one  week  to  10 
days  after  the  last  observation  used  in  the  orbit  fit.  These 
predictions  are  made  for  arbitrary  times,  and  the  calibration 
software  uses  numerical  interpolation  to  compute  the  observed 
quantities  at  the  time  of  observation.  The  interpolation  methods 
and  accuracy  assessment  is  given  by  Kantrowitz  (2).  It  has  been 
found  acceptable  to  use  a  time  spacing  of  1  minute  for  the 
predicted  data  points. 
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4.  Calibration  Model 

There  are  three  modes  of  tracking  with  the  Millstone  Hill 
Radar: 

1)  Coded  pulse,  used  primarily  for  acquisition. 

2)  Chirped  pulse,  with  up  to  1MHz  bandwidth,  with  the 
frequency  domain  tracker.  This  tracker  uses  coherent  integration 
to  improve  signal  to  noise  ratio,  and  a  fast  fourier  transform 

( FFT)  to  estimate  the  doppler  shift  for  estimation  of  line  of 
sight  velocity.  The  chirped  receive  pulse  is  deramped  in  analog 
and  sampled  with  the  matched  pulse  being  reconstructed  in  either 
a  hardwired  correlator  or  in  an  Array  Processor.  In  current 
operation  the  two  methods  should  generate  the  same  result. 

However  we  observe  the  correlator  and  AP  to  have  different  range 
calibration  constants. 

3)  Chirped  pulse,  with  up  to  1MHz  bandwidth,  with  the 
time  domain  tracker.  This  tracker,  designed  by  Raup  (1) 
provides  metric  data  on  non-coherent  targets  with  an  accuracy 
approaching  that  which  the  frequency  domain  tracker  achieves  on 
coherent  targets.  When  fully  implemented  one  can  expect  that,  by 
invoking  coherency  assumptions  where  appropriate,  significantly 
improved  range  rate  estimates  will  be  possible.  Since  the  time 
domain  and  frequency  domain  trackers  are  processing  the  same 
information  they  should  give  the  same  results.  That  means  any 
difference  must  be  software  induced.  Since  they  should  both  be 


9 


unbiased  estimates  of  the  same  quantity,  the  result  should  be  the 
same  to  within  the  formal  uncertainty.  This  does  not  seem  to  be 
the  case . 


Each  tracking  mode  should  have  its  own  mode  dependent 
calibration  model.  In  fact,  only  for  the  frequency  domain  data 
(2)  is  sufficient  tracking  data  in  hand  to  quantify  the 
calibration  constants.  Preliminary  data  is  available  for  the 
time  domain  tracker,  and  the  coded  pulse  and  correlator  data 
remain  to  be  investigated.  In  all  cases  a  basic  general 
calibration  model  is  adopted  with  some  mode  dependent 
parameters.  Therefore,  the  main  results  available  are  for  the 
chirped,  1MHz,  frequency  domain  tracker.  Preliminary  results  on 
the  other  trackers  (modes  of  operation)  will  be  given  when 
available.  The  basic  calibration  model  is  described  in  Table  II. 


Observable 


Range 

Range  Rate 
Elevation 

Azimuth 


Refraction 


TABLE  II 

CALIBRATION  MODEL 

Calibration 

Parameter 

Bias 

Bias-0 

Bias 

Droop 

Bias 

Traverse 

Misalignment 

Tilt 


Depends  On 


elevation 

elevation 
elevation , 
azimuth 
level  sensors 
elevation, 
elevation  rate, 
ionosond, 
ref ractometer 


Using  the  frequency  domain  tracker  we  have  found  that  range 
rate  data  is  unbiased  with  an  accuracy  of  1  mm/sec.  Table  III  is 
a  sample  of  three  weeks  of  range  and  range  rate  biases  for  June 
1984.  The  date  is  given  as  Modified  Julian  Date  (MJD)  which  is 
the  Julian  Date  minus  2400000.5,  and  day  of  year  ( DDD ) ,  and  the 
number  of  observations  is  n.  The  bias  for  each  pass  of  tracking 
data  is  determined  using  a  Lageos  reference  orbit  derived  from 
laser  ranging  data.  Figure  2  is  the  residuals  from  each  pass  of 
data.  Azimuth  and  elevation  residuals  are  also  given  for 
completeness.  The  unbiased  accuracy  of  1  mm/sec  (a  part  in  10 ') 
is  achieved  with  careful  treatment  of  the  tracking  system, 
refraction  and  the  theory  of  Special  Relativity  (4).  The  fact  of 
zero  range  rate  bis  is  so  generally  true  that  we  normally  do  not 
determine  range  rate  bias,  and  the  range  rate  bias  value  in  the 
calibratin  model  is  kept  at  zero.  Range  rate  is  therefore  the 
most  powerful  observable  in  determining  orbits  using  only  the 
Millstone  Radar  data.  With  an  accuracy  of  1  mm/sec,  the  range 
rate  is  measured  to  1  x  10" '  which  is  comparable  with  a  laser 
range  measurement  of  10  cm  (10"B). 

The  situation  with  azimuth  and  elevation  is  similar,  but  not 
so  promising.  Basically,  the  angle  measurements  are  limited  by 
mechanical  devices  (encoders)  to  0°.002.  This  is  a  respectable 
4x10" b  measurement.  The  calibration  model  for  the  azimuth  and 


TABLE  III 


MH 

CALIBRATION  WITH 

LASER 

REFERENCE 

ORBIT 

MJD 

DDD 

<R> 

n 

0 

<£> 

n 

a 

(m) 

(m) 

(mm/sec) 

(mm/sec ) 

45854.9 

155 

-4.6 

70 

3.0 

-1 .5 

71 

1.7 

45855.9 

156 

-7.3 

96 

1.5 

-0.6 

96 

3.1 

45857.9 

158 

-9.9 

25 

21.3* 

1.7 

32 

1.2 

45859.3 

160 

-13.5 

76 

4.1 

0.2 

75 

2.1 

45861.3 

162 

-10.5 

70 

1.7 

-0.2 

69 

2.1 

458863.3 

164 

-8.9 

83 

23.8* 

0.0 

88 

1.3 

45864.8 

165 

-8.6 

10 

1.8 

-0.5 

9 

2.1 

45865.2 

166 

-9.8 

98 

3.0 

1.0 

98 

2.6 

45866.3 

167 

3.7 

104 

5.3 

1.2 

104 

3.0 

45867.0 

168 

-8.9 

87 

3.1 

-2.0 

88 

2.9 

45868.4 

169 

-7.8 

60 

3.0 

-0.0 

61 

2.8 

45870.3 

171 

-3.7 

36 

2.3 

-1.7 

34 

1.5 

45872.3 

173 

-13.2 

84 

3.4 

1.1 

81 

3.8 

45873.2 

174 

-4.9 

83 

1.8 

1.8 

72 

2.6 

*  Coded  Pulse  Data 


TABLE  IV 


AZIMUTH,  ELEVATION  CALIBRATION  MODEL 


Calibration  Parameter 


Dec  1983 


February  1984 


Azimuth 

bias  (mdeg)  -45.0 

collimation  (mdeg)  -12.5 

Axis-skew  (mdeg)  -3.20 

Elevation 

bias  (mdeg)  290.0 

slope  (mdeg/deg)  -0.323 


-45.0 

-16.0 

-0.333 


293.0 

-0.347 


elevation  seems  stable  in  the  sense  that  over  periods  of  many 
months,  the  parameters  have  small  changes.  These  changes  appear 
to  be  seasonal,  but  too  little  data  is  in  hand  to  be  sure.  As  an 
example,  the  calibration  model  adopted  in  December  1983  and  in 
February  1984  are  given  in  Table  IV. 

These  azimuth  and  elevation  calibration  parameters  are 
directly  applied  to  the  data  going  into  the  metric  data  base,  and 
therefore,  no  further  correction  need  to  be  applied.  These 
parameters  are  given  here  for  illustration. 

There  is  no  evidence  that  the  model  parameters  for  azimuth 
and  elevation  are  different  in  the  three  tracking  modes  and  one 
would  not  expect  any  difference.  The  situation  is  made  less 
clear  because  the  angle  measurements  have  not  been  completely 
reliable,  i.e.  a  significant  number  are  discarded  in  any  orbital 
test.  This  provides  less  reliable  data  for  establishing 
calibration  parameters,  and  hence  less  certain  calibration 
values.  Further  discussion  of  angle  data  is  left  to  another 
report . 

The  range  calibration  for  each  pass  is  adequately  modelled 
as  a  bias.  Table  III  is  a  sample  for  three  weeks  (June  1984)  of 
the  bias  for  each  pass  determined  using  a  Lageos  reference  orbit 
determined  with  laser  ranging  data.  There  is  apparently  real 
pass  to  pass  variations  in  this  calibration  parameter. 

Discussion  of  these  data  will  be  given  in  section  6.  Suffice  it 
to  say  that  the  bias  determination  is  known  to  approximately  1 
meter:  2xl0-8  for  a  synchronous  satellite  or  lxl0“7  for  Lageos. 


5.  Refraction 


The  basic  measurement  is  a  two  way  time  delay.  The 
radiation  travels  at  the  speed  of  light  in  vacuum,  and  more 
slowly  in  the  atmosphere  or  ionosphere,  so  called  refraction.  To 
obtain  an  estimate  of  the  true  range  or  range  rate,  we  must 
estimate  the  contribution  of  refraction  on  the  observed  time 
delay.  The  refraction  can  conveniently  be  given  as  a  correction 
to  range  and  range  rate. 

5.1  Tropospheric  Refraction 

Many  analyses  of  refraction  modeling  exist,  and  the  basic 
relations  need  not  be  described  (5).  We  therefore,  describe  the 
observed  quantities  available. 

The  troposphere  extends  to  about  30km  altitude,  and  for 
radio  frequencies  less  than  the  22.5  GHz  water  vapor  line,  is 
nondispersive;  i.e.  the  index  of  refraction  (nt)  is  independent 
of  frequency.  Now  it  is  conventional  to  write 

nfc  =  1  +  Nfc  x  10- b 

where  N^  is  the  tropospheric  ref ractivity .  The  factor  10“b  is 
conventionally  included  and  results  in  Nt  values  of 
approximately  300.  For  the  troposphere,  N<-  depends  on  the 
composition,  pressure,  temperature,  and  relative  humidity  of 
air.  The  most  commonly  used  expression  for  an  average 
atmospheric  composition  is  (5): 
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where  P  is  the  pressure  in  millibars,  T  is  the  temperature  in 
degrees  Kelvin,  and  e  is  the  partial  pressure  of  water  vapor  in 
millibars.  Since  temperature  and  pressure  are  not  available  at 
altitude,  surface  data  are  used  for  the  surface  refractivity 
(N°t)  and  a  model  with  altitude  is  assumed: 

Nt(h)  =  N«  exp(-h/H) 

where  H  is  called  the  scale  height.  Average  values  of  N°t  and 
H  for  the  United  States  are  approximately  313  and  7km 
respectively.  From  the  on-site  ref ractometer  we  monitor  N°t. 
Therefore,  calculation  of  tropospheric  refraction  assumes  the 
scale  height  (H)  and  a  plane  stratified  atmosphere. 

5.2  Ionospheric  Refraction 

The  ionosphere  consists  of  ionized  particles  extending  from 
about  80km  to  beyond  1000km.  There  is  evidence  of  charged 
particles  (the  plasmasphere )  to  arbitrary  heights,  e.g.  to 
20000km,  (6).  We  must  assume,  for  the  moment,  that  these  higher 
charged  particles  have  a  small  effect  on  the  total  refraction, 
and  ignore  them. 

The  index  of  refraction  (ni)  is  less  than  1,  and  is 


expressed  in  terms  of  the  ionospheric  refractivity  as: 


nj  -  1  -  /  Neq2  \  -  1  -(^)2=  1-N. 

\e0Mw2  / 

where  q  is  the  electron  charge  1.6021x10“ 19  coulomb/electron 
M  is  the  electron  mass,  9.1091xl0“31  kg/electron 

ec  is  the  permitivity  of  free  space 

1/(4»x10“7c2) -10“9/( 36n)  farad/ra 
Ne  is  the  electron  density  (electrons/m3 ) 

f  is  the  frequency  of  the  propagating  wave  ( 2 w £— co )  (  Hz) 

fQ  is  equal  to  8.978/Ne  is  called  the  plasma 
frequency 

c  is  the  speed  of  light  c  =  299792458  m/sec. 

Since  n^  depends  of  frequency  (u),  the  ionosphere  is 
dispersive.  Now,  if  w<w0  n*  cannot  be  real  and  no  wave  will 
propagate,  i.e.  there  will  be  total  reflection.  For  typical 
values  of  Ne,  a>o  *  1°  MHz.  Therefore  the  L-Band  (1295MHz)  and 
UHF  (400MHz)  frequencies  are  well  above  the  plasma  frequency. 

With  an  ionosonde  or  backscatter  radar,  one  transmits  a 
number  of  frequencies.  Those  frequencies  below  w0,  are  all 
reflected.  Measuring  the  delay  of  the  reflected  wave  at  any 
frequency,  gives  an  estimate  of  the  height  of  the  electron 
concentration  corresponding  to  that  resonant  frequency.  The 
maximum  frequency  reflected,  therefore  corresponds  to  the  maximum 
electron  density  (Nmax)  and  is  known  as  fop2*  Now,  for 
calculation  of  refraction  in  range  and  doppler  measurements,  one 


30 


needs  the  integrated  electron  content  or  Total  Electron  Content 

(TEC)  along  the  propagation  path.  We  have  no  way  of  knowing  this 

from  the  bottom  side  ionosonde  measurements  since  they  tell  us 

nothing  about  the  distribution  of  electron  density  above  hmax, 

the  height  of  the  maximum  electron  density.  Neither  does  it 

indicate  anything  about  the  thickness  of  the  ionosphere.  It  is 

possible  to  use  a  general  model  of  the  ionosphere  and  employ  the 

observed  maximum  electron  density  to  scale  the  model.  It  is 

found  (7)  that  the  slab  thickness  S  =  TEC/N  is  rather  stable 

max 

and  well  behaved.  It  is  more  easily  modeled  than  either  TEC  or 

N _  ,  and  we  use  the  relation  that 

max 

TEC (instantaneous) =N  (measured) *S (model) . 

InaX 

Therefore,  based  on  these  arguments  (7)  we  use  a  slab  model 


which  is  equivalent  to  modeling  the  index  of  refraction  as 
(l  h  <  h_.„ 


N?  ={1  -  (foFj/f) 


h  <  h  <  h  +S 
max  max 

h  <  h  +  S 
max 


where  S  is  the  slab  thickness,  foF2  is  given  in  MHz  and  N  in 

3 

electrons  m  . 
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5.3  Refraction  Corrections 

The  index  of  refraction  (n)  allows  one  to  calculate  the 
transmission  time  of  electromagnetic  waves  from  the  propagation 
velocity  (v),  noting  that 
v  =  c/n 

where  c  is  the  speed  of  light  in  vacuo.  Strictly,  this  is.  the 
phase  velocity  (Vp).  For  a  nondispersive  medium  like  the 
troposphere  vp  <  c,  and  vp  is  the  same  as  the  group  velocity 
(Vg)  and  the  signal  velocity  (vs).  We  have  the  time  of 
transit  from  location  to  location  x2  as: 


Xl  Xx 

For  a  dispersive  medium,  such  as  the  ionosphere,  the  result 
depends  on  the  type  of  measurement.  In  the  ionosphere  the  phase 
velocity  vp=c/ni>c.  However,  the  group  velocity 
Vg=vp-dn/dw<c .  Note  that  for  this  dependence  of  n^  on  ai  we 
have  (vg  vp)=c2.  Now  for  measurements  involving  the  group 
velocity  we  have: 

x2  x2 

Tg  i  ^  *  5  /  (1  +  (1/2)Ni>  dX 

Xl  X! 


and  for  measurement  involving  the  phase  velocity  we  would  have: 


V 


^  -  i  j  »-« 


-(1/2)N.)  dx. 


Now,  one  cannot  measure  Tp,  but  it  is  possible  to  measure 
phase  differences.  Conceptually  the  difference  in  phase,  A$, 
between  the  return  signal  and  the  local  oscillator,  measured  in 
successive  pulses  separated  by  At  can  give  the  doppler  shift: 


Afi  = 


A*.-  A+i_1 


Now  using  approximations  sufficient  for  elevations  greater  than 
10  degrees,  the  additional  delay  due  to  refraction,  expressed  in 
meters  for  the  troposphere  is: 

Arfc  =  7  x  10“ 3  Nt/sin  E 

where  E  is  the  elevation  angle,  and  for  the  ionosphere,  using 
the  group  velocity,  is: 


“i  ■  k  p^y 


S  sec  x 


where  x  is  defined  by 

Ro 


sec  x  = 


'R2  -  R2  cos2  E 
o  e 


is  the  elevation 


Re  is  the  geocentric  station  radius 
R0=Re+hmax  is  the  geocentric  height  of  the 
maximum  electron  density  and 
S  is  the  slab  thickness. 


The  values  of  fQF2  are  measured  with  a  local  ionosonde, 
and  hmax  and  S  are  computed  from  models  by  (7).  These  models 
are  functions  of  solar  activity  (sun  spot  number  in  terms  of  the 


10.7  cm  solar  flux),  time  of  day  and  season.  In  case  Nt  or 
fQF2  are  n°t  observed,  seasonal  estimates  of  them  are  used. 
The  correction  to  range  rate  for  the  troposphere  is: 


-7x10" 3  Nfccos  E 
sin  ZE 


dE 

3t 


and  for  the  ionosphere,  using  the  phase  velocity  is: 


i  f  0^2  2  R*  2 

Ar  =  —  ( - }  S  sec  3  x  ( — )  cos  E  sin  E  — 

i  2  f  RQ  dt 


The  range-rate  measurement  from  the  frequency  domain  tracker  at 
Millstone  Hill  is  based  on  the  difference  in  phase,  and  hence  the 
correction  depends  on  the  phase  velocity.  For  the  time  domain 
tracker,  it  depends  on  the  group  velocity,  and  the  correction  has 
the  opposite  sign. 

As  an  example,  the  size  of  the  refraction  corrections  are 
estimated  as  follows:  Typical  values  of  the  troposphere  are 
N(-=313.  Therefore  we  have 

at  zenith  E=90°  Art  =  2. 2  m 

and  at  E=20 0  Art  =  6.4  m 
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For  the  ionosphere  with 


*0*2 

f 

S 

^max 


Ro 


12  MHz 

1200  MHz  (L-band) 
280  km 
300  km 

6678  km 


we  have : 


Elevation 

X 

Ari 

(°) 

(°) 

(m) 

90 

0. 

14. 

20 

63.8 

32  . 

and  with  the  same  assumptions  as  before,  we  have  for  the 
ionosphere  at  E  =  20° 


Ar^  =  1.5  cm/ sec. 

An  example  of  ionospheric  and  tropospheric  refractions  is  given 
in  Fig.  3. 

The  troposphere  correction  should  be  accurate  to  a  few 
percent.  The  ionosphere  correction  is  based  on  a  framework  of 
assumptions  and  approximations.  It  is  clear  that  the  slab  model 
is  a  significant  oversimplification  of  the  actual  ionosphere.  I 
is  expected  to  be  no  better  than  25%.  During  periods  of 
quiescent  ionosphere,  i.e.  away  from  day  night  transition  and 
solar  disturbance,  this  model  should  be  better  than  10%.  During 


6.  Orbital  Results 

The  basic  orbit  computation  program  is  a  numerical 
integration,  with  variable  step  size,  automatic  truncation  and 
round  off  error  control.  The  force  model  used  is  described  in 
Table  V,  and  the  reference  system  is  described  in  Table  VI.  The 
orbital  fit  to  laser  ranging  data  is  not  yet  at  the  accuracy  of 
the  data;  however  for  the  radar  data  the  orbit  accuracy  does 
exceed  the  data  accuracy. 

For  analysis  of  Millstone  Radar  Data  alone,  the  parameters 
described  in  Tables  V  and  VI  can  be  used  directly  for  orbit 
determination.  The  error  committed  in  this  way  is  less  than  the 
observation  uncertainty.  However  for  the  precision  laser  range 
data,  this  is  not  true.  In  particular  the  following  quantities 
must  be  determined  with  each  trajectory  (typically  8  days 
duration ) ; 

1 )  Drag  Parameter 

2)  Earth  Polar  Motion 

3)  Earth  rate  of  rotation. 

In  addition  there  may  be  laser  data  from  new  stations  whose 
coordinates  are  not  so  well  known.  In  this  case,  the  solution  is 
improved  by  including  the  three  station  coordinates  in  the 
solution  vector.  Therefore,  typical  laser  data  runs  have  as 
solution  parameters:  6  orbital  constants  +  drag  acceleration  +  3 
earth  rotation  parameters  +  3  station  coordinates. 


Table  V 


Force  Model 

Geopotential 
sun  &  moon 
Body  Tides 

Solar  Radiation  Pressure 
Earth  Albedo 
Drag 

General  Relativity 

Table  VI 
Reference  Frame 


GEML2  to  degree  and  order  20  (9) 
USNO  Lunar  Ephermeris  (10) 
k2  *  0.285  e  =  0.0° 

K  ■  1.002 
as  0.32 

p0  =  3.43x10“  ia  grams/cm 3 

Y  ■  1.0 


Precession  &  Nutation  IAU  J2000  System  (10) 

Polar  Motion  and  UTl  BIH  Rapid  service,  USNO 

Predictions  and  Determined  with 
Laser  Data 

Laser  Site  coordinates  (11) 


Station  Coordinates 
Tide  displacement 


h  =  0.60 


When  sufficient  data  from  a  new  station  is  available ,  reliable 
coordinates  can  be  determined ,  and  that  site  can  be  held  fixed. 
Figure  4  shows  the  range  residuals  for  a  pass  from  Millstone/  and 
the  range  residuals  from  an  overlapping  pass  of  laser  ranging 
data.  This  is  a  good,  but  typical,  result  from  an  8  day  orbit 
fit.  On  Fig.  4  are  listed  additional  solution  parameters  and 
their  formal  uncertainty. 

The  drag  parameter,  polar  motion  and  earth  rotation  data  are 
real,  time  dependent,  geophysical  information,  uniquely  provided 
by  laser  ranging  to  Lageos.  These  geophysical  data  will  be 
discussed  elsewhere. 

The  accuracy  of  the  Lageos  orbit  fit  using  quick  look  laser 
ranging  data  is  at  present  50  cm  r.m.s.  This  is  due  to  lack  of 
final  calibration  and  editing  of  the  laser  data,  and  the 
preliminary  coordinates  used  for  some  of  the  stations.  The 
accuracy  of  the  Lageos  orbit  fit  using  Millstone  radar  data  is  2 
to  5  meters  r.m.s.  The  larger  error  occurs  when  the  amount  of 
data  decreases.  For  extrapolation,  (modes  2  and  4)  the  errors 
grow  to  1  m  and  10  m  for  the  laser  derived  and  radar  derived 
orbits  respectively. 

The  determination  of  Range  Bias  has  been  made  using  mode  1 
since  December  1983.  The  values  are  given  in  Appendix  A  and  are 
plotted  in  Fig.  5.  These  values  are  corrections  to  the  nominal 
calibration  model  used  for  data  in  the  Metric  Data  Base. 
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Therefore  these  values  can  be  used  to  correct  data  in  the  metric 


data  base.  This  correction  can  be  applied  to  any  satellite 
observed  during  the  same  tracking  session  with  the  same  tracking 
mode.  For  each  entry  the  date  is  given  in  three  ways: 

1)  Modified  Julian  Day  (MJD  *  Julian  Date  -  2400000*5)  to  0.1 
days,  2)  calendar  date,  and  3)  day  of  year.  The  range  correction 
is  given  in  meters.  The  number  of  orbit  fits  containing  the  pass 
of  data  and  an  uncertainty  are  also  given.  When  the  number  of 
estimates  is  two,  the  value  given  is  the  mean,  and  the 
uncertainty  is  half  the  difference.  For  more  than  two  estimates, 
the  mean  is  given  and  the  uncertainty  is  the  root  mean  square 
difference  from  the  mean.  Also  given,  when  available,  is  the 
range  bias  derived  from  a  reference  orbit  computed  from  laser 
ranging  data.  A  recent  history  of  the  radar  operation  is  evident 
from  this  figure,  and  the  annotations.  For  example,  the  change 
in  calibration  accompanying  installation  of  a  new  transmitter 
component  is  dramatic.  The  rapid  detection  and  consequent  change 
of  the  nominal  calibration  model  was  possible  because  of  the 
routine  use  of  Mode  2  calculation  of  range  bias. 

These  calibration  data  are  strictly  for  the  frequency  domain 
1  MHz  chirped  tracker.  The  calibration  correction  given  in 
appendix  A  can  be  added  to  data  in  the  metric  data  base  for  any 
satellite  taken  with  the  frequecy  domain  tracker  during  the  same 
tracking  session.  Range  data  corrected  in  this  way  can  be 


considered  to  have  an  accuracy  associated  with  the  bias  value 
used  i.e.  2  to  5  m.  This  has  been  done  with  success  for  a  number 
of  investigations. 

There  are  some  tracks  taken  with  the  experimental  1  MHz  time 
domain  tracker.  The  range  bias  values  are  consistent,  but 
significantly  different  from  the  others  in  the  time  period. 
Apparently  the  time  domain  tracker  is  not  consistent  with  the 
frequency  domain  tracker  from  the  range  calibration  stand  point. 
More  data  are  required.  Clearly  it  is  desirable  for  the  two 
trackers  to  have  the  same  calibration;  the  frequency  domain 
tracker  calibration  result  can  then  be  applied  to  the  time  domain 
tracked  data,  and  only  one  calibration  model  need  be  maintained. 

It  must  be  emphasized  that  this  calibration  result  does  not 
apply  to  the  so  called  coded  pulse  data.  The  coded  pulse  data  is 
a  narrow  band  alternative  tracking  mode.  The  signal  to  noise 
ratio  for  Lageos  results  in  range  precision  of  80  to  120  meters. 
At  present  it  is  used  during  the  acquisition  phase  of  all  passes 
and  a  number  of  points  observed  this  way  are  in  the  metric  data 
base.  It  can  also  be  true  that  a  whole  pass  is  taken  in  coded 
pulse  mode.  However  so  little  data  is  available  on  Lageos,  that 
nothing  can  be  said  about  the  calibration  of  the  coded  pulse 
data. 

A  small  number  of  tracks  were  taken  with  the  correlator 
mode  rather  than  the  array  processor  mode.  There  is  strong 


evidence  of  roughly,  a  15  m  bias  introduced  using  the 
correlator.  To  demonstrate  this  a  track  was  taken,  where  the 
tracking  mode  was  switched  in  mid  pass  from  correlator  to 
frequency  domain,  and  then  to  time  domain  tracking.  The  range 
residuals  are  plotted  in  Fig.  6,  which  illustrates  two  points. 
First,  the  frequency  domain  tracker  with  the  correlator  has  a  15 
meter  bias,  with  respect  to  the  frequency  domain  tracker  and  the 
time  domain  and  frequency  domain  trackers  have  appear  to  a 
different  range  calibration  bias. 


Residuals  for  correlator,  frequency  domain,  and  time  domain  tracking. 


7.  Prediction  Accuracy 

The  results  reported  here  are  essentially  from  an 
posterori  analysis.  To  move  to  Mode  2  or  Mode  4  operation,  which 
would  provide  more  immediate  calibration  values,  requires 
predictions  of  the  satellite  position  and  the  earth's  position. 

As  mentioned  above,  precision  orbit  determination  requires 
simultaneous  determination  of  the  satellite  state  vector  and  a 
number  of  time  variable  geophysical  parameters.  The  geophysical 
parameters  that  most  stress  the  system  are  the  "along  track"  drag 
like  force  and  the  UT1.  The  variability  of  both  of  these  is 
apparently  due  to  meteorological  effects  and  are  not  easily 
predictable.  We  give  in  Fig.  7  the  along  track  acceleration 
recovered  from  the  laser  ranging  data.  In  Fig.  8  is  given  the 
run  off  of  predicted  UTl,  provided  by  the  U.  S.  Naval  Observatory, 
from  the  actual  UTl  values,  provided  by  the  Bureau  International 
de  l'Heure  (BIH) .  In  both  cases,  most  of  the  time,  the  prediction 
error  is  within  5  meters  after  10  days.  Investigation  of  the 
best  methods  of  predicting  these  phenomena  is  just  beginning. 

At  this  stage,  these  provide  the  theoretical  limit  to  use  of  Mode 
2  or  Mode  4  for  calibration.  This  limitation  is  true  for  any 


satellite  based  calibration  system. 
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8.  Summary  &  Conclusions 

A  method  for  Radar  Calibration  has  been  described  and 
demonstrated.  It  is  apparent  that  routine  1  m  a  posteriori 
calibration  is  now  available  when  suitable  observing  proceedures 
are  used.  It  is  anticipated  that  this  will  be  available  in  real 
time. 

There  are  three  important  elements  in  the  observing 
procedures 

1)  the  tracking  mode  must  be  known, 

2)  sufficient  number  and  distribution  in  elevation  of  data 
points  must  be  available,  and 

3)  refraction  measurements  must  be  made. 

A  posteriori  calibration  can  be  improved,  and  for  the  most 
precise  data  analysis,  these  calibration  corrections  can  be 
applied.  All  the  basic  modes  of  tracking,  (coded  pulse, 
frequency  domain,  time  domain)  with  and  without  the  correlator 
have  different  range  calibration  models.  Azimuth  and  elevation 
have  the  same  model  for  each  mode.  The  frequency  domain  range 
rate  is  an  unbiased  observable  with  an  accuracy  of  1  mn/sec.  The 
bias  in  the  time  domain  range  rate  has  yet  to  be  assessed. 

Since  September  15,  the  nominal  system  calibration  has  been 


within  ±  5  m,  for  the  frequency  domain  tracker.  One  can 
speculate  that  many  of  the  earlier  outliers  were  due  to 
unrecorded  system  configuration  changes.  Also,  the  general 


variability  of  system  calibration  has  decreased  with  time,  and 
possibly  procedures  are  more  faithfully  followed.  Certainly 
system  calibration  variability  may  be  due  to  a  personal  equation. 
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